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Abstract 

We discuss the properties of two Bose-Einstein condensates in different spin states, represented 
quantum mechanically by a double Fock state. Individual measurements of the spins of the particles 
are performed in transverse directions (perpendicular to the spin quantization axis), giving access to 
the relative phase of the two macroscopically occupied states. Before the first spin measurement, the 
phase is completely undetermined; after a few measurements, a more and more precise knowledge of 
its value emerges under the effect of the quantum measurement process. This naturally leads to the 
usual notion of a quasi-classical phase (Anderson phase) and to an interesting transposition of the 
EPR (Einstein-Podolsky-Rosen) argument to macroscopic physical quantities. The purpose of this 
article is to discuss this transposition, as well as situations where the notion of a quasi-classical phase 
is no longer sufficient to account for the quantum results, and where significant violations of Bell type 
inequalities are predicted. 

Quantum mechanically, the problem can be treated exactly: the probability for all sequences of 
results can be expressed in the form of a double integral, depending on all parameters that define the 
experiment (number of particles, number and angles of measurements). We discuss the differences 
between this case and the usual two-spin case. We then discuss the effect of the many parameters 
that the experimenters can adjust for their measurements, starting with a discussion of the effect of 
the angles of measurement (the "settings"), and then envisaging various choices of the functions that 
are used to obtain violation of BCIfSH inequalities. We then discuss how the "sample bias loophole" 
(often also called "efficiency loophole") can be closed in this case, by introducing a preliminary 
sequence of measurements to localize the particles into "measurement boxes". We finally show how 
that the same non-local effects can be observed with distinguishable spins. 

1 Introduction 

Two of Einstein's many famous contributions to physics are the theoretical discovery in 1925 of the 
Bose-Einstein condensation in a gas ,1^ and, ten years later, the celebrated Einstein-Podolsky-Rosen 
(EPR) argument 2J. Both have had an enormous influence in the discipline and stimulated much work, 
both theoretical and experimental. Although they are both fundamental, these contributions appear 
almost completely disconnected: the former is more "standard" physics, initiating the domain of quantum 
statistical physics, with many applications in gas and condensed matter physics; the latter belongs to 
the foundations of quantum mechanics and has indeed attracted the attention of philosophers. It is 
therefore interesting to realize that both contributions are connected logically, and that the study of 
interfering Bose-Einstein condensates may shed a new light on the fundamental debate initiated by EPR. 
The basic reason is that, while the original EPR argument applies to two microscopic particles, its 
transposition to Bose-Einstein condensates introduces systems that can be macroscopic. The "elements 
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of reality" , introduced by EPR as attached to microscopic particles, then characterize the relative phase 
between the two condensates which, in the case of spin condensates, may determine a macroscopic spin 
orientation. This is an important difference: while one can argue as Bohr that elements of reality 
for microscopic particles do not exist independently of the measurement apparatuses, it is more difficult 
to deny that macroscopic physical systems possess an independent physical reality. As a consequence, 
the EPR argument then becomes more compelling j^. Another interesting feature is that, while in the 
original example the EPR element of reality (and therefore an additional, or "hidden" variable) appears 
as completely foreign to standard quantum mechanics, with condensates the relative phase emerges rather 
naturally within its formalism of standard quantum mechanics, simply as a consequence of particle number 
conservation [4]. 

A natural question then is whether condensates can lead to violations of local realism, as systems 
of two particles do, in other words whether violations of Bell-type inequalities [51 [B] also occur with 
macroscopic condensates. In view of this macroscopic character, one could expect that the answer to this 
question is no. Actually, it turns out that it is yes: as shown in a letter T recently, strong violations do 
occur when measuring the individual transverse spin orientation of particles in a double spin condensate 
with equal populations, even for arbitrarily large condensates. An essential condition, nevertheless, is 
that the spins of all particles be measured; missing one or more spins makes the violation disappear. 
The relative phase that emerges under the effect of the first few quantum measurements behaves like a 
quasi-classical variable, so that violations are impossible in this regime; nevertheless, if one continues the 
sequence of measurements until its maximum, one reaches situations that can no longer be understood 
with this classical phase, but recover an intrinsic quantum character. 

The purpose of the present article is to discuss in more detail several aspects of the questions treated 
in [7] as well as those that were not be treated there. In §[21 we define the physical system under consid- 
eration and recall the notation, as well as previous results. In §[31 we come back to the EPR argument in 
the context of the relative phase of two condensates, and in § [H to violations of Bell inequalities obtained 
within stochastic local realist theories (as opposed to the usual two-spin case where deterministic local 
realist theories are more natural). Generally speaking, a useful feature of transverse spin measurements 
in condensates is that one can calculate exactly the effect of many parameters defining the measurements: 
angles at which the spins are measured, number of particles and number of measurements, various func- 
tions used to obtain violations of the inequalities, etc.; the discussion of the effect of these parameters is 
given in §[5] with numerical analysis. The next section, §[^1 deals with the well-known "sample bias loop- 
hole" (often also called "detection, or efficiency, loophole" ) , which can be closed (in thought experiments) 
with the help of preliminary measurements, exactly in the perspective proposed by Bell [8J to handle such 
situations. In § [3 we show that the violations of local realism are not related to boson statistics, but 
that the same violations can be obtained with distinguishable spins, provided they are in an appropriate 
initial state. We then draw conclusions in § [SI 

2 Summary of previous results 

We recall in this section the results already obtained in [4J [T] , starting with the equations that show how 
a phase naturally emerges from the predictions of standard quantum mechanics in a series of transverse 
spin measurements. We then discuss in what conditions this phase behaves as a quasi-classical quantity, 
or retains a strong quantum character; in the latter case, we briefiy introduce the violation of Bell 
inequalities and local realism that can be obtained, keeping a more detailed discussion for §§[3] and [HI 

2.1 Quantum predictions 

Consider an ensemble of 7V+ particles in a state characterized by an orbital state u(r) and a spin state 
and iV_ particles in the same orbital state v{r) with spin orientation — . This physical system is described 
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quantum mechanically by a double Fock state: 



$ > 



(«",+) 



vac. > 



(1) 



where au.+ and a„._ are the destruction operators associated with the two populated single-particle states 
and |vac. > is the vacuum state. The total number of particles is: 



N = N+ + N^ 



(2) 



The operators associated with the local density and the local density of spins can be expressed as function 
of the two fields operators ^'±(r) associated with the two internal states ± as: 



n(r) = *^^_(r)4'+(r) + 5'^_(r)*_(r) 
a,(r)- *^.(r)*+(r)-*L(r)*_(r) 



^'!_(r)^'+(r) - *^(r)^'_(r) 



(3) 



while the spin component in the direction of plane xOy making an angle ip with Ox is: 

a^(r) = e-''^'J'^r)*_(r) + e''^*^r)*+(r) (4) 

When a measurement of this component performed at point r provides the result i] = ±1, the corre- 
sponding projector is: 



Pri=±i{r,(p) = ^ [n{r) + Tj a^{r)] 



(5) 



For a series of N measurements that are performed at different points (ensuring that the projectors 
all commute) along directions ipi, the probability for obtaining a series of results r/i ± 1 can be written as: 



V{r]i,r]2,...r]N) = < * I (^l, </'l) X ^'»J2 (i"2, (^2) x •••-P'?iv , <^Ar) I $ > 



(6) 



Strictly speaking, this expression is not a probability but a density of probability, which must be integrated 
in a finite volume Ar to provide a probability; for instance, the probability for finding a particle in volume 
Ar is: 



P(Ar) 



n{r ) 



*^(r')1'+(r') + ^l{r')^^{r'] 



(7) 



A similar expressions for finding its spin along any direction </> is obtained by integrating ^ over space. 
Here, for convenience we do not write the integrals explicitly, but it must be understood that expression 
([6]) as well as those we write below are in fact integrated over all position variables ri, r2, etc. in spatially 
disconnectecG volumes Ai, A2, etc. As in the second ref , we call Ai, A2, .. the "detection boxes" 
and refer the reader to this reference for a more detailed discussion, including the case where no particle 
at all is detected in the detection box; we also come back to this point in § [S] 

As in ref. [7] we now substitute the expression for aip{r) into (O and 0. In the product of projectors 
appearing in ([H]), all r's are different and commutation allows us to push all the field operators to 
the right, all their conjugates to the left; then that each ^'+(r) acting on | $ > can be replaced by 
u(r) X au,+ , each 4'-(r) by v{r) x au^~, and similarly for the Hermitian conjugates. With our initial 
state, a non-zero result can be obtained only if exactly iV+ operators au,+ appear in the term considered. 



^If there was an overlap between the "detection boxes", the projectors would no longer commute and equation (|6j would 
no longer be valid; it would then be necessary to use the complete "Wigner formula" with twice as many projectors ordered 
in a symmetric way. 
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and -/V„ operators a^,^-; a similar condition exists for the Hermitian conjugate operators. To express 
these conditions, we introduce two additional variables. The first variable A ensures an equal number of 
creation and destruction operators in the internal states ± ; each \['_|_ (or au,+) is multiplied by e'^, and 
each vj/^ (or aj, _(_) by e^'"^ (operators related to the — spin state remain unchanged), and the conservation 
of the number of particles in spin + states is then obtained by the following integration: 



inX _ r 
T;- e — Onfl 



(8) 



(since the total number of particles is unchanged, the number of — particles is then also automatically con- 
served) . The second variable A expresses that the difference between the number of destruction operators 
in states + and — is exactly 7V+ — -/V_ , through the integral: 



7;— e e ^ + ' — On^N 



(9) 



This time, each (or a„,+) is multiplied by e~ and each (or a^^_) by e (but creation operators 
remain unchanged). The introduction of all these exponentials into the product of projectors ([5|) in 
then provides the expression (c.c. means complex conjugate): 



N 



N+m^l n 9 ^^'^ + \vir,)f + V, U*(r,>(r,)e'(^-'^^+'^) + c.c 



(10) 



where the field operators are replaced by a constant number Af-|_!iV_! because, after integration over A 
and A, the only surviving terms are all associated with the same matrix element, that of the product of 



iV-|_ operators aj, ^ and iV_ operators a|j _ followed by the same sequence of destruction operators inside 
state I $ >. To simplify the equations, from now on we assume that all measurements are made at points 
where the two wave functions u(r) and u(r) are equaj^; we can then write the probability as: 



c.c. 



(11) 

By changing one integration variable (A' = A + A) and using A parity, we obtain: 

1 f^'^ dA f^^ dX' ^ 

P(77i,772,...w) =^jv^ y — COS [{N+-N^) A] J _ J] [cos (A) + 77, cos (A' - v'j)] (12) 

where the normalization coefficient Cn is obtained by writing that the sum of probabilities of all possible 
sequences of ry's is 1 (we come back to this point in §[6]): 

— cos[(7V+_iV_)A] [cos(A)]^ (13) 

In the above equations, it is sometimes convenient to reduce the integration domain of A to the interval 
between — 7r/2 and +Tr/2. This is possible, since (before variable A is changed into A') changing A into 
vr - A introduces into ^ a factor („i)W+-^-+^^ or (-1)^^+, which is 1. 



If they have equal modulus but different phases, their phase difference at point rj is simply added to the measurement 
angle ipj. 
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Let us assume unequal populations, for instance Nj^. > N^. Then, in the product over j contained 
in (jl2p . only some terms provide a non-zero contribution in the integral over A; in fact, we must choose 
at least N+ — N- factors contributing through cos (A), and therefore at most N — (iV+ — A^_) = 2N- 
factors contributing through the rjj and (fij dependent terms. We see that 2N^ is the maximum number 
of spins that can provide a transverse spin measurement that depends on the result 77 and on the angle 
of measurement; all the others have equal probabilities 1/2, whatever the angle of measurement is. This 
is physically understandable, since (A^-i- — N^) spins + are unmatched with — spins, and can thus not 
be found in a coherent superposition of the two spin states; they then provide 1/2 probabilities for any 
direction of transverse spin measurement. We therefore see that all spins can contribute coherently to 
the measurement only if — iV_. 

The above equations are valid only when the number of measurements M is equal to N , but can 
still be used when it is smaller. The reason is that any sequence of M < measurements can always 
be completed by additional A^ — M measurements, leading to probability (HH). We can therefore take 
the sum of over the results of the missing N — M measurements, which amounts to setting the 
corresponding rjj equal to zero and doubling the remaining term. We then obtain: 

1 ^iA /'+7r 

V{m,V2, -mi) = J ^ cos [(A^+ - N_)A] [cos A]^-"' J — [] [cos (A) + ry, cos (A - ^,)] 

(14) 

Finally, what we will need below is the value of the quantum average of the product of results, i.e., the 
sum 7^(771,772, .•.77m) times this product over all possible values of the 77's. This sum, according to p^ . is 
given by: 

E{^,,^2,-Vm) ^{Cn)-' —cos[{N+~N.)A][coaAf-'' ^J|cos(A-^,) (15) 



2.2 Classical or quantum regime of the phase 

The relative phase A is a priori completely absent from Fock states in standard quantum mechanics, 
but nevertheless appears naturally in its formalism as a mere consequence of the conservation of the 
number of particlef|j. The phase actually occurs in an integral; for the first measurement, the integral 
merely expresses that the phase is initially completely undetermined, as one could expect; for a series 
of measurements, the phase integral provides the relation between the successive results and introduces 
their correlations. Depending on the measurements, this phase takes a classical or a quantum character. 

If the number of experiments M is much less than a very large A^, and if A^+ — because cos A^~*^ 
peaks up sharply at A = 00, equation (fH)) becomes: 

"'-^ j=l 

We then recover the results of [4] as well as of previous work (refs [9] to [E]). For a given A, the proba- 
bilities can be obtained by considering that the sample is completely polarized in a transverse direction 
determined by angle A; the spin measurements then become independent processes with individual prob- 
abilities given by {1 -I- 77^ cos (A — tpj)}, exactly as for a single isolated spin; the additional ingredient is 
the A integral, which expresses that all values of A are equally probable and introduces the correlations. 

^More precisely, the conservation of the difference of the number of particles in the two internal states, which is the 
conjugate variable of the relative phase. 

''Here we take the point of view where the A integration domain is between — 7r/2 and +7r/2; otherwise, we should also 
take into account a peak around A = tt. 
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In this case, all predictions of quantum mechanics lead to predictions that are perfectly compatible with 
the idea of a pre-existing phase, which takes a well-defined value before any measurement and remains 
constanllfl; this value is initially completely unknown and changes from one realization of the experiment 
to the next. This fits well with the concept of the Anderson phase, which originates from spontaneous 
symmetry breaking at the phase transition (Bose-Einstein condensation): at this transition point, the 
quantum system chooses a phase, which takes a completely random value, but then plays the role of a 
classical variable in the limit of very large systems. 

On the other hand, if TV — M is not a large number, the peaking effect of cos A^~^^ does not occur 
anymore and A can take values close to 7r/2, so that the terms in the product inside the integral are no 
longer necessarily positive; an interpretation in terms of probabilities then becomes impossible - see § 
14.21 for more details. In these cases, the phase does not behave as a semi-classical variable, but retains a 
strong quantum character; the variable A controls the amount of quantum effects. 

2.3 Violations of local realism in the quantum regime 

Indeed, in the quantum regime where A is not limited to values around zero, equations (|14p and (jlSp 
may contain strong violations of Bell inequalities and therefore of local realism. Consider a thought 
experiment with two condensates, each in a different spin state (two eigenstates of the spin component 
along the quantization axis Oz). The two condensates extend into two remote regions of space Da and 
Db where they overlap and have equal orbital wave functions, and where two experimentalists Alice and 
Bob measure the spins of the particles in arbitrary transverse directions (any direction perpendicular to 
Oz), see Fig.[TJ We assume that all measurements performed by Alice are made along the same direction 
ifa, which plays here the usual role of the "setting" a, while all measurements performed by Bob are 
made with another single angle iph- To complete the analogy with the usual situation with two particles, 
we assume that Alice retains of all her measurements just their product A, while Bob retains only the 
product B of his results (these numbers are both ±1, the parities resulting of all the local measurements) 
- other possibilities than these simple products are considered in § [5l 

We now assume that Alice, in successive realizations of the experiment, uses two possible orientations 
If a a.nd if'g^, and that Bob does the same with two possible orientations iff, and ip'f^. Within local realism, 
for each realization of the experiment, it is possible to define two numbers A, A', both equal to ±1, 
and associated with the two possible products of results that Alice will observe, depending of her choice 
of orientation; the same is obviously true for Bob, introducing B and B' . Since AB ~ ±A'B', either 
AB + AB' or A'B — A'B' vanishes, and the following inequality is straightforward: 



As a consequence, the average of the products in (|17p . obtained by repeating the experiment and the 
measurements many times, must be between —2 and +2; this is the famous BCHSH (Bell, Clauser, Horne, 
Shimony and Holt) inequality [121 HO] j a consequence of local realism. 

Within standard quantum mechanics, the above reasoning no longer holds. As emphasized by Peres 
[H], "unperformed experiments have no results", so that several of the numbers appearing in (fT7|) are 
undefineqj; in fact, only two of them are defined after the experiment has been performed with a given 
choice of the orientations. Consequently, while one can calculate from (llSp the quantum average (Q) of 
the sum of products of results appearing in the middle of ([T7)) . there is no special reason why (Q) should 
be limited between +2 and —2. Situations where the inequality: 



^In our calculations, we have assumed that the measurements arc all made in a very short period of time, so that the 
evolution of the system between them can be ignored. 

more precise statement would be; "for a given realization of the experiments, unperformed experiments have no 
results; or, if they do, each result depends not only on the local setting but also (and non-locally) on the other remote 
setting" . 



2<AB + AB'± {A'B - A'B') < 2 



(17) 



2< (Q) <2 



(18) 
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Figure 1: Two condensates with fixed number of particles, one having up spins and the other down spins, 
overlap in two remote regions of space Da and Db where measurements of transverse spin orientations 
are made by Alice and Bob. In the upper part of the figure, we assume that each of the orbital wave 
functions associated with the two condensates are simply connected, and extend continuously between 
the two remote measurement sites. This is not a necessary condition, nevertheless. For instance, as shown 
schematically in the lower part of the figure, we can assume that the wave function of each condensate is 
coherently split into two parts and separates into two disconnected components (each of which will then 
contain a fluctuating number of particles) , and that each of which propagates to one measurement region 
separately. 

is violated are sometimes called "quantum non-local situations" . 

Situation in which (|16p holds can not lead to such violations, since this equation contains positive 
probabilities inside the integral and has precisely the form from which Bell inequalities can be derived (we 
come back to this point in more detail in 14. ip . So, here we no longer assume that M <^ N but consider 
the other extreme, M = N . The simplest case occurs when iV+ = = 1 and when Alice and Bob make 
one measurement each; it is then easy to see that ([l} defines a triplet spin state | 5 = 1, Ms = >. In 
this case, it is well-knowr0 that a violation of occurs, by a factor \/2 when the angles form a "fan"[f| 
spaced by x = 7i'/4; this saturates the Cirel'son bound [22]. But the violations also occur for arbitrarily 
large values of and 7V_: for instance, consider the case M — N, Na = 1 (Alice makes one measurement 
only) and = N — I (Bob makes the maximum number of remaining possible measurements). Then 

''When one measures the components of the spins along directions that are perpendicular to the quantization axis, the 
predictions of quantum mechanics are the same for this triplet case and the singlet state | 5 = 0, Ms = >, provided one 
just reverses the direction of one measurement. 

*The term "fan" refers to the angles arranged as ipab = Vba' = Vb'a = X ^-nd v'6'a' = 3x where ip^b = Va ^ Vb- 
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equation ([15]) becomes: 





1 -1 



dX' 



(19) 



— TT 



— TT 



27r 
2tt 



1 



cos((Pq - V^b) 



This result is precisely the same as that for the two-particle case, so that the Cirel'son bound is saturated 
again. For other values of Na and Nt,, substantial violations of continue to occur [7]; we refer to §[5] 
for more details. 

3 A macroscopic element of reality 

We begin this section with a brief discussion of the Leggett-Sols argument which, from a different point of 
view, leads to the same conclusions as the EPR argument; we then discuss in more detail the transposition 
of the EPR argument to spin condensates. 

3.1 The Leggett-Sols argument 

Leggett and Sols ^23j discuss a situation that has strong similarities with double spin condensates: two 
superconductors, initially in Fock (number) states, are coupled by a Josephson junction, which creates a 
current flow between them; the phase of this time oscillating current corresponds to the relative phase of 
the two superconductors. In standard quantum mechanics, initially this phase is completely undetermined 
but, as soon as the time dependence of the current is measured, the phase is created by the very act of 
measurement; it takes some random, but well-defined, value. The authors ask "Does the act of looking to 
see whether a Josephson current flows force the system into an eigenstate of current, and hence of relative 
phase?" . The answer of standard quantum mechanics to this question is "yes" , but the authors point 
out that "if one thinks about it seriously, this answer is bizarre in the extreme" . To illustrate why, they 
suppose that the current is of order of, say, kiloamps, and that it is measured through the observation 
of a small magnet needle. They then ask: "Can it really be that by placing a minuscule compass needle 
next to the system, with a weak light beam to read off its position, we can force the system to realize 
a definite macroscopic value of the current? Common sense certainly rebels against this conclusion, and 
we believe that in this case common sense is right" . In other words, because the current is arbitrarily 
large, its phase can not be created by a tiny measurement apparatus; it must already have existed before 
the measurement. Since this "element of reality" is not contained in standard quantum mechanics, this 
theory is incomplete. Here the argument is not locality, as in the EPR argument, but simply that a very 
small system can not completely modify an arbitrarily large system through a quasi instantaneous and 
mysterious measurement process, without any precise physical mechanism to explain why and how. 

Double condensates undergoing transverse spin measurements, with a number of measurements M <gC 
N, are very similar to the case discussed by Leggett and Sols. From (fTC|) . we can obtain that the probability 
of finding r]M in the M-th measurement along ipk , having found r?i , • • • r/M-i in the previous measurements 
along (fii^ ■ ■ -ipM-i- 





where: 



9mW ^ Y[ {l + ??j Cos(A- v?j)} 
i=i 



(21) 
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(these equations are valid only if M ^ N). The evolution of this "probability function" was studied in 
ref. [25] where it was shown that it peaks up sharply after only a few measurements. A typical result 
is shown in Fig. [2] where the measurements are done at a single angle Lpa = 0. Two peaks arise because, 
with a single angle of measurement, the sign of the relative phase is not determined, but a very small 
number of additional measurements at a new angle causes the collapse to a single peak. 




Figure 2: The angular distribution g{^) as a function of angle for three different numbers of measurements 
of transverse spin, 10 measurements (dashed line), 150 measurements (dottted line), and 300 measure- 
ments (solid line). For a single measuring angle this always has two equal peaks, corresponding to the 
ambiguity of the spin direction with respect to the transverse plane; but making measurements along 
another direction rapidly removes one of the peaks. The peaks narrow when the number of measurements 
increases. 

The transposition of the question of Leggett and Sols would then be "Can it really be that, by 
measuring the transverse direction of a few microscopic spins, we can force the macroscopic polarization 
of 10^^ atoms (or more) to take a definite value?" The analysis of double spin condensates in the present 
paper shows that, within standard quantum mechanics, one can obtain detailed and exact predictions 
of the effect of an arbitrary number of measurements in any direction on the macroscopic polarization; 
the paradox can then be studied in more detail [4]. The subject is also related to the discussion of 
spontaneous symmetry breaking of spin condensates given by Siggia and Riickenstein [24] . 

3.2 Transposition of the EPR argument to double condensates 

One usually discusses the standard EPR argument in the form proposed by Bohm, with two spin 1/2 
particles entangled in the singlet spin state 15 = 0, Ms = 0), or equivalently in the triplet |S' = 1, Ms — 0). 
In these situations, the measurement of the spin of the first particle in any direction determines the value 
of the second spin along the same direction. As soon as the first result is known, the second is also known 
with certainty when the two directions of measurement are parallel: perfect correlations are predicted by 
quantum mechanics. This leads EPR to their famous statement : "If, without in any way disturbing a 
system, we can predict with certainty the value of a physical quantity, then there exists an element of 
physical reality corresponding to this physical quantity." 

With condensates, what emerges from the measurements (still assuming M ^ N) is a relative phase 
(Anderson phase) of the condensates, through the process discussed above. Initially, this phase is com- 
pletely undetermined, and the first spin measurement provides a completely random result. But the 



9 



phase rapidly emerges under the effect of a few measurements, and then remains constanlQ; it takes a 
different value for each realization of the experiment, as if the experiment was revealing the pre-existing 
value of a classical quantity. 

Assume now, as in |3], that the double condensate extends over a very large region of space, covering 
both Alice's laboratory and Bob's very remote laboratory, as shown in Fig. [TJ We then have a situation 
where, without in any way disturbing Bob's system, we can predict from Alice's results the direction 
of the macroscopic orientation that Bob will observe in his remote laboratory; then there must exist in 
Bob's laboratory an element of reality associated with this prediction - at this stage, standard quantum 
mechanics still agrees, provided one uses the postulate of wave packet reduction, which accounts for 
this element of reality. In addition, since Bob's laboratory is far away and therefore protected from 
any influence of Alice's operations, the element of reality also necessarily existed before Alice made any 
measurement - then standard quantum mechanics cannot agree anymore. In fact, it does not only ignore 
this initial element of reality, but even says that the phase is completely undetermined before the first 
measurement; the EPR argument then concludes that quantum mechanics is incomplete. We have already 
emphasized in the introduction that the major difference between this case and the usual two spin case 
is that, here, the EPR element of reality can be macroscopic; this weakens Bohr's rebuttal of the EPR 
argument, which hinges on the ambiguity of physical reality for isolated microscopic systems (considered 
independently of the macroscopic measurement apparatuses), and seems more difficult to transpose to 
macroscopic systems. 

There are also a few other differences. First, with double condensates, a single measurement of the 
spin of one particle is not sufficient to determine the relative phase; Alice and Bob, have to measure the 
spin of at least a few particles to obtain a reasonable determination of this phase, with better and better 
accuracy when the number of measurements increases. This is not a problem, since the total number 
of available particles may be macroscopic, while a few tens of measurements are already sufficient to 
obtain an excellent determination; see [21] for a discussion of the strategies that Alice and Bob may 
use to optimize their knowledge of the phase. We remark in passing that, with condensates, the usual 
discussion of incompatible measurements, counterfactuality, etc. is not relevant: Alice and Bob can 
use exactly the same experimental procedures in all realizations of the experiment, and obtain a good 
knowledge of the phase. 

The second difference is that, while for two particles the quantization axis along which both spins 
polarize is fixed by the direction of first measurement, here the system "chooses for itself its phase and 
therefore its quantization axis; the emerging transverse orientation can have any direction with respect 
to the direction of measurements. Moreover, this direction is only known with an accuracy that is limited 
by a quantum uncertainty, which decreases when the number of measurements increases (phase/number 
quantum uncertainty relation). Even if Alice and Bob choose parallel directions for their measurements 
(or any relative direction), perfect correlations are not predicted in general by quantum mechanics, but 
only equal probabilities for obtaining result +1 for instance; individual measurements therefore remain 
stochastic processes so that, strictly speaking, the words "with certainty" used by EPR therefore do not 
apply with double condensates. 

Fortunately, this does not ruin the EPR reasoning: if Alice makes appropriate measurements to 
determine the phase with good accuracy, and if Bob chooses a direction of measurement that is parallel 
to the spontaneous transverse magnetization that Alice has observed, the certainty is just replaced by a 
high probability, 99% for instance. Alternatively, one can also consider that Alice and Bob use sequences of 
individual spin measurements to measure the angle of the transverse spin polarization; if these sequences 
are sufficiently long, there is a high probability that their determinations of the phase will agree within 
a small error bar. Therefore strong correlations are still obtained in this case, even if no longer at the 
level of individual measurements. Local realism then ascribes their origin to correlated elements of reality 
belonging to these remote regions of space, and the essence of the EPR reasoning still applies. We note, 

®As already mentioned, we ignore any evolution of the system between measurements. 
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nevertheless, that here only one additional element of reality emerges from the EPR reasoning, that 
associated with the direction that the system has chosen, while in the usual situation with two spins 
all components of Bob's single spin are predictable from Alice's result (provided she chooses a parallel 
direction of measurement). But, if one accepts local realism, one missing element of reality is already 
sufficient to prove that quantum mechanics is incomplete ! One can summarize all this discussion by 
saying that the usual EPR microscopic elements of reality, associated with all components of a single 
spin, collapse here into one single, macroscopic, element of reality. 

4 Microscopic violations of local realism 

We now continue the EPR reasoning to derive Bell inequalities; we then show that the quantum predic- 
tions violate these inequalities; we complete this section with a comparison between violations obtained 
with GHZ and double Fock states. 

4.1 Bell inequalities within stochastic local realist theories 

The derivation of Bell inequalities from the EPR conclusions involves different reasonings in the usual 
case (two spins) and for two condensates. We first recall the situation in the usual case. 

4.1.1 Two spins 

With two spins, the derivation of the Bell theorem starts from the existence of well defined functions 
A{X, ifa) and B{\, ipt) giving the results of the measurements; these functions depend on the fluctuating 
elements of reality A that each particle carries with it, and of the local orientation if a or ipb of the 
measurement apparatus. Within local realism, their existence is proved by the fact that, for any direction 
chosen by Alice (or Bob), it is always possible that Bob (or Alice) will choose a parallel direction; one can 
then predict with certainty the second observed result from another measurement made very far away. 
The results of spin measurements are therefore deterministic functions of the additional variable A and 
of the local setting; the original Bell reasoning O E] then leads to the usual Bell inequalities. 

The result can be generalized to a stochastic point of view; the inequalities do not require determinism, 
but can also be proved within stochastic reahst theories, provided they are local [20]. We call -P+(A, (fa) 
the probability that Alice will obtain a result +1 when the relative phase is A and when she has chosen a 
direction ipa for her measurements, P^{X,Lpa) the probability for the opposite result; a similar notation 
P±{X, <fb) is used for Bob. For a given reahzation of the experiment, with a given phase A, the expectation 
of the product of the results is: 

Pl{\, Va)Pi{\ ^b) + P- (A, (^a)P- (A, ^b) - Pl{\ Va)P- (A, ^fc) - P1{X, <fa)Pl(X, fb) (22) 

The average of the product of the results observed by Alice and Bob in many realizations of the experiment 
is then: 

— [PU\Va)-Pl{X,^a)][PUX,^b)-PH\^b)]; (23) 

that is, the average over the possible values of the relative phase A of the product of the two quantities: 

A(A, ^a) - (A, Va) - P^(A, ^a) = 2P^{X, ^a) - 1 

B{X,^b) ^ PliK'Pb) - P^{\^b) = 2PliX,^b) - I ^ ' 

In the second equation of each line, we have taken into account that the sum of probabilities P+ and P_, 
for given A and angle of measurement, is 1; because all probabilities are numbers between and 1, for 
any value of A and the angle Lp both A and B are numbers between —1 and +1. 
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If now we form the combination of the average that appear in (|17p . we obtain the average over the 
phase A of the expression: 

AB +AW±(A'^~A'B') (25) 



where the primes indicate that the angle ipa has been replaced by (p'^^ (or ipi, by ip'f^). The only difference 
with the deterministic case is that the numbers that appear in ((25|) are no longer equal to ±1, but 
have some value between —1 and +1. But expression (|25p is linear with respect to all of these numbers 
separately. Therefore, if we replace one of the numbers, A for instance, by ±1, we change the expression 
to new values that provide upper and lower bounds of the initial value. Doing the same thing for all four 
variables in succession therefore provides new lower and upper bounds which, since now all numbers are 
±1, are ±2. At the end of the process, we see that (|25|) is still bound between —2 and +2; its average 
value over the phase A must have the same property, so that the BCHSH equations remain valid. 

4.1.2 Two condensates 

With two condensates, the situation is different: when and N- are more than 1, equation 
does not contain situations with full correlations. When individual spin results observed by one of the 
experimenters cannot be predicted with certainty (for any direction of measurement) from the result 
already obtained by the other, local realism can no longer be used to derive the existence of functions 
A{X, ipa) and B{X, (pb). There is no way to force the axis of quantization, as already discussed in ? 13.11 the 
relative phase that emerges from the measurements is independent of the directions of measurements. If, 
for instance, Alice and Bob choose a common direction that happens to be perpendicular to the transverse 
direction that has spontaneously appeared, each of them will have 50% probabilities for the two results, 
and the EPR element of reality provides them with no information at all. The connection between the 
EPR reasoning and the Bell theorem can therefore not be directly transposed from the two-spin case. 

A way to proceed is to extend the analysis of § 12.21 by a reasoning that we will call the quasi-classical 
treatment of the relative phase (Anderson phase). We have seen that the local realist EPR argument, 
applied to sequence of measurements where M remains smaller than the number of particles N, leads 
us to conclude that the sample is fully polarized in some unknown direction. This full polarization has 
no reason to disappear when more measurements are performed: for instance, if Bob's sample is initially 
fully polarized, it will keep this full polarization if Alice accumulates more measurements on her side, 
and even completes the sequence of measurements so that M becomes equal to N: arbitrarily remote 
experiments can not influence the local physical properties of Bob's sample. So we can consider, within 
local realism, that both Alice and Bob actually do experiments on fully polarized samples with unknown 
transverse directions. In this case, for each realization of the experiment, all spins are in the same 
individual quantum state, and the spin measurements are actually independent processes. We can then 
write the simple formula: 



and where the correlations between the measurements are introduced by the A integral in (|26p . Of course 
the simplest idea is to choose for all of them the same probability, for instance that given by quantum 
mechanics for the measurement on a single spin: 




(26) 




(27) 



- [1-f ?7 cos (A - V3)] 



(28) 
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but we can also take for this probability an arbitrary function of its variables, provided conditions (j27p 
are fulfilled. In any case, we arrive at a situation where the EPR reasoning leads to probabilities instead 
of certainties. 

Assume now that Alice makes Na measurements and Bob Ni,; formula gives the probability of 
any series of results they observe. For using the BCHSH formula, both must choose functions A and B 
that depend on their local results, and take values that remain between ±1. There is a large flexibility at 
this stage: Alice decides to attribute value +1 to some some chains of her results ryi, 772, ■• VNa, value 
^ = — 1 to all the others; Bob makes a similar choice. Now, to obtain the probability that the product 
AB is 1, we can sum the probabilities of two exclusive events (either A = B = 1, or^ = i3 = —1) and 
use (f26| . Let us for instance calculate the probability of the first event which, according to ([26|l , is the 
product of two local probabilities, P°(A;(/3i, --(pNa) a-nd P^{X; (pN^+i, ..ipN^+N^), defined as: 

P^(A;<^i,..(pArJ= ^ n^-if (^''^J-) ; Pl{X;m.+i,-^N.+N,)^ n ^^i'(A,^,) (29) 

A=+lj=l B=+l j=Na + l 

where the two sums are taken over the sequences of ry's that realize A — I for Alice, S = 1 for Bob; 
similar reasonings apply to the other values of A and B. At this point, we see that we have made the 
connection with the previous calculation: we can define functions A and B by replacing in (|24p A hy A 
and B by B, and the rest of the reasoning goes unchanged; the only difference is that each local angle 
ipa and ipb is replaced by a series of angles. The essential property remains: each function A or B still 
depends only on the local angles chosen by its experimenter, and the BCHSH inequality is still valid. 

We conclude that the quasi-classical treatment of the relative phase leads to the BCHSH inequalities; 
each time we can write the probability of combined measurements in the form (|26|). where the P^^ {(pj) 
are numbers between and 1, these inequalities hold. 



4.2 Comparison with the quantum predictions 

The quantum predictions of equation (fH|) are not exactly of the form (|26p . but they are similar. Are the 
differences sufficient to introduce violations of local realism? We already know that they are, since we 
have seen in S I2.3l that the BCHSH inequalities can be violated by the quantum results; here we study in 
more detail the mechanism of these violations. 

First, we have already noted that, if N is large and if M ^ A^, the peaking effect [cosA]^~^^ selects 
only the values of A around zero, so that a good approximation is to take A = inside all the brackets 
contained in the product over j; then the A integral disappears and one exactly recovers (|16p. so that no 
violation is possible. For large violations of local realism, the most interesting cases occur when M has 
its maximal value A'^; so, while in the preceding section we discussed mostly the situation where M <^ N, 
here we are mostly interested in the opposite case. 

If A does not vanish, pi)) remains similar to (pS)) . while not identical. The first difference is that 
contains a double integral, but this is not essential: clearly the results of S I4.1l can easily be generalized to 
more than one additional variable A, for instance two A and A, and to situations where the distribution 
p(A, A) is not uniform; any positive normalized distribution is possible. If we attempt to bring p4p to a 
form that is compatible with local realist theories, we must satisfy conditions (|27p : for this purpose, in 
the product over j, we factorize [cos A]*^ so that each term in the product becomes: 



cos (A — ip-j) 
■' cos A 



(30) 



Then the sum of probabilities associated with the two rj ~ ±1 results is indeed 1, as requested. But, at the 
same time, we see that the "probabilities" introduced in this way may become negative for some values of 
the variables, which opens the way to violations of the BCHSH inequalities, by a mechanism that we now 
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discuss. First, if we define A and B as equal to 1 for any value of the variables 77, because the quantum 
probabilities are normalized to 1 by summing over the results 77's, each of the 4 terms contained in < Q > 
is then exactly 1, so that < Q >— 2. Now suppose that, for some values of the variables, the product 
of "probabilities" is negative; if we redefine A and B in such a way that makes their product negative 
for these values, this will automatically increase the value of < Q > beyond 2 and violate the BCHSH 
inequalities. In fact, since A and B depend only on the ry's and not on the other variables A and A, this 
operation may affect at the same time domains of the variables where the product of "probabilities" is 
negative, and positive; the net effect is then a balance between positive and negative contributions, and 
the violations do occur when the contribution of the former outweighs those of the latter. In any case, 
negative probabilities are a necessary conditions for violations of the inequalities; in § [U we discuss in 
detail several examples of these situations. 

4.3 GHZ states versus double Fock states for violations of local realism 

Mermin pSj has proposed a thought experiment involving many particles and leading to exponential 
violations of local realism; we now briefiy compare his scheme with ours. He uses a maximally entangled 
spin state (GHZ, or NOON state), which is sometimes considered as the "most quantum state" accessible 
to an N particle system. The GHZ states are also sometimes called "Schrodinger cat states" , since they 
involve a coherent superposition of states that are macroscopically distinct if N is very large; they are 
not easy to produce experimentally with many particles - to our knowledge, the world record [27] for 
the number of particles is iV = 5 - and very sensitive to decoherence ^2B] . 

Our double Fock state ([T|) is the simplest possible state that is compatible with Bose statistics. 
Conceptually, there is no simpler way to put together identical particles in two different spin states; at 
first sight, it does not even look entangled but, still, strong violations of the BCHSH inequalities do 
occur. Reference |29j shows how double Fock states with equal populations can be used in interferometers 
to measure the relative quantum phase at the Heisenberg limit. Such states also undergo decoherence 
by coupling to the environment, although more slowly than GHZ states [3H]; the "natural basis" for 
decoherence is given by phase states (corresponding to different macroscopic spin orientations), and its 
effect on our conclusions are minor, since nothing in the calculations requires coherence between various 
phase states. With present experimental techniques, there seems to be no enormous difficulty in producing 
double Bose-Einstein condensates. Nevertheless, to observe the quantum non-local effect we study here, it 
is essential to obtain the equality — iV_ . This means for instance that it is necessary to carefully avoid 
atom losses in both condensates: values of and 7V_ of the order of 10 seem accessible experimentally, 
but probably not orders of magnitude more. 

A striking feature of Mermin's thought experiment is the exponential violation that is predicted; we 
obtain nothing similar here, just a violation comparable to Cirelson's limit. Nevertheless it should be 
realized that, for N particles, the observable that Mermin introduces is the sum of 2^"^ commuting 
products of operators. It seems difficult to imagine how to measure this sum without measuring the 2^~^ 
commuting components. Seen in this way, the Mermin scheme amounts to accumulating 2^~^ measure- 
ments, and taking a sum of results in a way that accumulates the violation and makes it proportional to 
the square root of the number of measurements; this procedure can of course be implemented in other 
schemes, including two-particle experiments, or our scheme, and leads to violations that are even linear 
in the number of measurements. But the price to pay in all cases is a large increase of the number of 
measurements. 

Finally - and this is probably the most important difference that we have already emphasized — there 
is an important conceptual difference, since with the GHZ state the EPR "elements of reality" remain 
microscopic, while with double Fock states they may be macroscopic. 



14 



5 Types of measurements; numerical results 



Here we consider various types of measurement and the quantitative values of violations that occur with 
Bose-Einstein condensates. While up to this point we have considered only the usual form of inequality 
shown in (fT7|l . other forms are possible as we see here. Moreover, the quantities A, B, etc. can take on 
forms other than a simple product of ry's. Values of the inequality violations will be given in this section 
but angles for the spin measurements will be presented in an Appendix. 



5.1 Products of ?7's 

We return to (jlSp in which we computed the average of a product of experimental results for the r]i. Such 
a product is ±1 and so qualifies to be an A or B. We consider the case where the number of experiments 
is equal to the total number of particles, M — N and the numbers of up and down spin particles are the 
same, iV_(- = . Then we have the simple result: 

LJf ncos(A-^,) 

We assume Alice makes P measurements, all at the same angle tpa, and Bob makes — P at angle (fib, 
corresponding to products of results A and B; A' and B' correspond to two other values of the angles 
and (^1,. Equation ([3T|) then reduces to: 



„, ^ J-TT 27r — V" -ray-"" V" T-oy 



f cos ^(A-y,)cos^ -^ (A - 
/riflcosA] 

Then the quantum average for the Bell test of the inequality p7p is: 

(g> - E{ipa, ipb) + E{ip'^,^b) + E{^a. V'b) ~ Eiip'^,ip',) (33) 
In the numerator of p2p we change variables to A' = A — iph; if we define: 

X = (pa-'Pb (34) 

we obtain: 

cos'^ (A' — x) cos^^^ (A') = / cos^^^ (A') [cos Acosx + sinAsinx] 

27r 27r 



fc=0 

The integral is known and we find: 



^, , (A^/2)!^^^ P!(A^-2A:)! . p-2k 

E[x) ^ ' > jrr sm^^YCOs^ y (36) 

N\ ^^^k\{P-2k)\{^-k)\ ^ ^ V ; 



where {P/2} is the integral part of P/2. This result is efficiently evaluated in numerical maximizations 
of {Q), so that rather large A^ values can be treated. We always find a fan arrangement with ^a — ^b = 
'^b — '^a' = fb' ~ Va = X ^^'^ '/'b' ^ '/'a' ~ 3x, although thc value of x a-t maximum decreases with 
increasing A^- see the Appendix. 
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For arbitrary N and P = 1, we again see that E{x) — cosx as noted above in p^ . yielding (Q)^ 
2V2. For P = 2, the resuh is: 



TV- 1 



+ 1 



N - 1 



cos 2x 



(37) 



so that (Q)rriax ^^^^ depend on N. For iV = 4 the value is 2.28, but for very large N we obtain (Q)„jax ~ 
2.414; surprisingly it increases with iV. In Fig. [3] we plot {Q)^^y. versus N for various P values; even in 
the case P = N/2, we still get violations for all N. 
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Figure 3: The maximum of the quantum average (Q) for Alice doing P experiments and Bob N ~ P, as 
a function of the total number of particles N. Local realist theories predict an upper limit of 2; large 
violations of this limit are obtained, even with macroscopic systems {N -^00). Not shown is the case 
P = 1 for which the Cirel'son limit of 2\/2 is obtained for all N. 



In the case in which both P and N ^ P are very large, a simple approximation for E is available. By 
expanding the logarithm of cos^ x to second order in x we obtain the approximate form: 



cos y = e 2 y 



from which wc find: 



J —00 



P(W-P) 2 
2N ^ 



r+°° dXe 

J —OQ 



(38) 



(39) 



In the case P — N/2 we maximize a sum of Gaussians, and again find a fan arrangement, with a maximum 
of 8/(3 X 3^/®) = 2.32 at x = ^In3/A^. So the fan opening decreases as 1/^/N; this approximate result 
is valid for N as small as 12. This is similar to the results obtained by Drummond [BIT, who nevertheless 
starts from a many particle quantum state that is very different from ours (he considers N identical pairs 
of particles, each pair in a state that involves 4 single-particle states) as well as a different measurement 
scheme (number of bosons in one spin state, at two different locations). 
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We have also relaxed the constraint that all the angles within a set of measurements providing A, 
or those providing etc. be identical. Thus we might introduce ipai,'Pa2, etc. and ipi,i,Lpb2, etc. We 
have numerically examined such generalizations for N up to 10, and found that the maximization in this 
larger space collapses to the one we discuss above, ipai — 'fa2 ^ ■ ■ ■ ^ fa, with a single angle for all A 
measurements, and a single one for all B measurements, etc. 

Experiments with M < A'^ do not result in violations of the Bell inequalities. To see this write equation 
in the form: 

M 

/-:£ncos(A-^,) 
E{^uf2,-V^M) = 7+^-^^ — G{M,N+,N^) (40) 

where, if M is even: 

G(M,iV+,iV_)- /+;|Aeos[(iV+-W_)A]cos«A . . 

_ (JV-M)!M!Ar+!Af„! V^"^/ 

- (iV+-f )!(7V_-f )!(f .)^iV! 

and G = if 7\f is odd. The first factor in equation ([^0]) is the expectation value for M experiments with 
M particles and can lead to a violation of the Bell inequality. But the correction term G can be shown, 
for fixed M, to be largest for 7V+ = = N/2; then an analysis of G{M,N/2,N/2) shows that this 
quantity is always less than or equal to 2/3 unless M = N. Since (2/3)2V2 < 2, we have the remarkable 
result that one must measure every particle's spin in order to see a violation of the Bell inequality. Even 
missing the measurement of one or two particles ruins the observations of the quantum effect. 

5.2 Other definitions of A and B 

In the above analysis, we have used only a product of all the rj's as the A or B quantity. Other possibilities 
are available. For example, if Alice and Bob each make N/2 measurements, we might take A and B in 
the form: 

V1+V2 + --- Vn/2 ^^2) 
\vi +f?2 H VN/2\ 

which, if both Alice and Bob choose one single angle of measurement, would be a macroscopic polarization 
of the spins measured by each; more precisely, the numerator of this expression is the macroscopic 
polarization (in dimensionless units), and the denominator ensures a "binning" operation that retains 
only the sign ±1. Averaging the product of the two polarizations for (AB), we found that this procedure 
does not lead to a Bell violation for any set of angles, except of course for = 2 where the violation is 
2^2. For TV = 4, we find (Q)^^^ = 1.88; for iV = 8, we find 1.78; for N = 10, we find 1.970; for N = 14, 
we find 1.966; the value seems to converge asymptotically to 2, that is, the upper limit of local realism. 

On the other hand, if Alice makes — 1 measurements and Bob just one, the average of the product 
of her polarization and his single value: 

7?1 +7?2 H ?7jV-l , , 

does lead to violations for one value of N only. The values for iV = 4, 6, 8, 19 are respectively, (Q)„jax ~ 
1.41, 2.121, 1.59, 1.99. There is a violation for N = 6 with higher values again possibly tending to 2. 
This violation is a remarkable result since, here, Alice makes a measurement that is almost mesoscopic; 
the quantum character of Bob's measurement is nevertheless sufficient to maintain a significant violation 
of local realism. Nevertheless, if the number of measurements made by Alice increases beyond = 6, the 
violations disappear. Other cases where Bob makes two or more measurements, and Alice the complement 
to A^, do not lead to violations. 
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We have also tried considering averages of sums of two measurements in the form of products of 
ivi + ^2)/2. Such a quantity has possible values 1,0,-1 but the Bell inequality still holds. We found no 
cases where the quantum average of such pair averages lead to a violation. 

5.3 Other inequalities 

Other inequalities besides that of (fT7|) are possible. For example, consider the inequality: 

- 2 < ^{AB + A'B + AB' - A'B')iCV + C'V + CD' - C'V) < 2 (44) 

where each of the letters represents an 77 or a product of any number of ry's. We assume that the angles 
of measurement corresponding to the 77's in each letter are all the same (but releasing this constraint 
in the corresponding quantum average does not increase the violation found.). For = 4 each letter 
represents just one ry, while for N ^ 8 each letter corresponds to the product of two 77's. For TV = 6, 
A would represents one 77 and B two, etc. We find the violations shown in Table I when maximizing 
the corresponding quantum averages (the N ~ 00 result in Table I comes from making a Gaussian 
approximation, as in equation psp . for the powers of cosines in the integrals). 

Table I. Results for the {Q)^^^^ corresponding to the inequality of (jU) 

N {Q) 
4 2.66 
6 2.33 
8 2.18 
12 2.17 
00 2.15 

An extension of the idea in ([H]) is the inequality: 
-2 < ^{AB + A'B + AB'-A'B'){CV + C'V + CV'-C'V'){£J' + £'F+£J''-£'J'') < 2 (45) 

The quantum counterpart of this yields 2.66 for A'^ = 6; we therefore obtain a large quantum violation of 
this particular inequality, which therefore provides an interesting generalization of the BCHSH inequality. 
We have been able to treat larger N values in this case only by use of the Gaussian approximation discussed 
earlier. The violation continues for larger N with a limit of (Q)i„ax = 2.09. 

6 Sample bias (efficiency) loophole 

Our quantum calculations are consistent only if the " measurement boxes" are spatially disconnected (to 
ensure commutation of the quantum field operators) and if their volume A is sufficiently small to limit 
the number of particles in each of them to or 1; otherwise, the expressions of the projectors we use are 
not valicF°l. This means that the average number of particles in each box is much less than one, so that 
most measurements detect no spin at all. But, if one counts for all these non-detection events, clearly 
the quantum average of the product of results becomes very close to zero, and no violation of the Bell 
inequalities remains possible! 

This is not an unusual situation in Bell-type experiments. When detecting pairs of photons for instance, 
most photons are lost because of the finite solid angle that is captured by the detectors and of their limited 

^"The expression we use is actually the total amount of spin orientation within the volume (in h/2 units). If two particles 
are found inside the same volume A, then the two possible values of this orientation are ±2 (instead of ±1 for one particle), 
leading to eigenvalues 2 or of the (space integrated) operator This is in contradiction with the eigenvalues 1 or of a 
projector; moreover, values exceeding 1 are in contradiction with the assumptions leading to the Bell inequalities. 
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quantum efficiency. This is known as the "sample bias loophole", "efficiency loophole", "pair selection 
loophole" , etc. To avoid the problem, what is done in practice by experimentalists is to redefine the sample 
of events in the calculation of the averages: instead of the sample of all emitted pairs, they consider the 
sample of pairs for which particles in coincidence are indeed detected. This restores the possibility of 
a violation of the inequalities, but at the same time destroys the validity of the BCHSH inequalities 
themselves, since local realism stricto sensu is then no longer sufficient to derive them. The reason is 
that there is no way to ensure that the sample remains independent of the settings of the apparatuses, 
while this assumption is crucial for the proof of the Bell theorem: if the "settings" introduce a bias in 
the sample, the distribution of variables A may depend on them, and the proof of the theorem becomes 
impossible One then has to introduce extra assumptions, for instance that the measured probability is 
the product of a probability of detection (independent of the settings) by a spin (and setting) dependent 
probability that is relevant to the Bell inequality violation. This experimental loophole has been pointed 
out many times, and some authors (a minority) have even refuted all locality experiments for this reason; 
there is a large amount of literature on the subject. 

Fortunately, at least for thought experiments, the loophole can be closed; it is therefore not a 
fundamental obstacle, but only contingent on our present technologies. John Bell had a elegant way to 
solve the problem [5] , with the introduction of either "veto detectors" or "spin independent preliminary 
detectors" ; the purpose of these detectors was to properly define a sample of systems that is independent 
of the settings and ensures that a spin signal is always obtained at each detector. Similarly, Clauser and 
Shimony |^D] introduce "event ready detectors" that have the same function. In our case with Bose- 
Einstein condensates, we need something similar. The simplest idea is to assume that, before any spin 
measurement, spin-independent detectors are used to ensure that one particle (and one exactly) is found 
in each measurement box. It may be necessary to repeat the preparation procedure many times before 
this result is obtained, since in most cases no particle is found in at least one of the measurement boxes, 
but in theory this is not a problem: it is sufficient to ignore these cases, and to repeat the preparation 
stage as many times as needed until the desired result is observed. Only after this sample preparation 
stage has been successful will the spin measurements be performed. Alternatively, one can decide to 
replace the initial quantum state, the double Fock state, by the new state obtained after wave packet 
reduction is applied after a positive preparation stage. Here we study the explicit form of this new initial 
quantum state. 



6.1 Calculating a new quantum state 

The initial double Fock state is given by ([T]); we now assume that M "measurement boxes" are defined 
in the volume occupied by the orbital wave function, and decompose this wave function as: 

M+l 

u{r) ^ ^ u,n{r) (46) 

m— 1 

where for m < M the function Um(r) is the normalized "projection"!^ of the wave function u(r) into the 
measurement box number m; iiM+i(r) is defined as the complementary projection of u(r) outside all the 
measurement boxes. The Xm are the components of the linear decomposition of u(r) onto the contents of 
the various boxes, with: 



AI+l 



(47) 



It is even possible to show that a selection of detected pairs that is dependent on (f>a and (pf, makes it possible to 
reproduce any correlation of results within local realist models. In other words, local realism does not introduce inequalities 
anymore. 

^^This projection is equal to u{r) within the box, zero outside, and then normalized to 1. 
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For m < M, the smaller the measurement boxes, the larger the values of the normalized u„i(r)'s inside 
their boxes, and the smaller the coefRcients Xm', for m = M + 1, xm+i then remains close to 1. The 
creation operators for the states | u, + > and | u, — > can be expressed as functions of the creation 



operators a^^ _,_ and a^^ _ for the the states 



> and 



>: 



m — 1 



M+l 



al 



(48) 



so that: 



M+l 
m— 1 



al 



AT, 



Af+1 

E 

m —1 



vac) 



(49) 



All operators in this expression commute, so that usual algebraic expansions of the powers of sums can 
be used without special care. 

We see in (|49p that the double condensate state vector contains components on states where the 
number of particles in each box varies substantially. If, for instance, one selects inside both sums only 
the terms m = M + 1, all particles go to the complementary box, while all measurement boxes remain 
empty; no particle can be detected at all by any of the apparatuses. If, on the other hand, one selects 
only terms corresponding to a given measurement box, all particles accumulate into this particular box, 
while all the others remain empty. Of course, one can also spread the particles among all measurement 
boxes, and what we wish is to consider situations where they are equally filled. 

We then decide to introduce a new initial state by retaining from (|49p only the components where each 
measurement box contains one particle exactly. In other words, we project | $ > onto the subspace where 
each measurement box contains one particle, and obtain a new state vector | $ >. Mathematically, this 
vector could be written with the introduction of exponentials and integrations into (j49p , but for simplicity 
we do not write this expression; what is important here is not so much the mathematical form of the new 
initial state | $ > , but the fact that it exists and can be built, without changing the relative probabilities 
calculated with | $ > for detecting single spins in each measurement box. This is true by construction: 
one can easily see that that all components of | $ > that have been eliminated from | (f> > play no 
role whatsoever in the calculation of the probability for single particle detection; they just eliminate 
no-detection and multiple detection events. The only difference is the normalization of the ket, which is 
changed by the removal of all these useless components; the remaining components must be increased to 
restore normalization. Physically, since we are now sure that one particle, and one particle only, will be 
detected in each box, we know that the probabilities for all possible results rjm = ±1 add to 1, which is 
exactly what we wish for a violation of the inequalities. 



6.2 Discussion 

Expression ((49|) is nothing but the product of A''-|_ sums associated with internal state + by sums 
associated with the other with internal state — . One can see the physical system as the juxtaposition 
of two entangled subsystems, one corresponding to the content of the measurement boxes (system I), 
the other to the content of the complementary box (system II). But, if we limit the discussion to the 
case where M has its maximal value iV+ + iV_ (we have seen in § [5] that this corresponds to maximal 
violations of the BCHSH inequalities), then system II becomes empty; its quantum state is independent 
of that of the M measurement boxes and factorizes out. System I is then in a pure state, with components 
that depend on how the two internal states are distributed among the measurement boxes; the number 
of possibilities is: 

iV! _ A^! 

{N -N+)\N+\ ^ {N-N^)\NA ^ ' 
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In state | > , by construction, both Alice and Bob perform local experiments on a fixed total number 
of particles, but with a fluctuating number of spins + and spins — . This remark allows us to come 
back to the situation shown in the lower part of figure [T] When the wave function of each condensate 
is coherently split into two disconnected parts, each spin system remains a single condensate with an 
orbital function that is the coherent sum of two distant components, and the number of particles in each 
of these components has large fluctuations. Therefore, neither Alice nor Bob knows the number of spin 
up, or spin down, that she/he receives; the corresponding fluctuations are essential for the interesting 
quantum non-local effects to occur. The total number of particles contained in Alice's sample also has 
large fluctuations in state | $ >, but since this state can be replaced by | $ > without affecting our 
results, we see that these fluctuations are not essential; in state | $ >, the fluctuations between the 
numbers of spins up and spins down in each region of space are correlated in a way that cancels the 
fluctuations of the total number of particles in this region, without destroying the non-local effects. 

Rather than writing the resulting state for system I in more detail with identical particles, it is simpler 
to consider now distinguishable particles. 



7 Violations of local realism with distinguishable particles 

We now show that our results for violations of the Bell inequalities are not just limited to Bose-Einstein 
condensates but also apply to distinguishable particles. When the number of measurements M has its 
maximal value N = + A^_, state | $ > corresponds to all particles localized in different boxes, with 
no spatial overlap, so that they do behave as distinguishable particles; if we wish, we can number them 
in the state vector by assigning them the number of the box they occupy; this operation does not affect 
the physical predictions. 



7.1 Quantum state 

Using the numbering of the boxes in which the particles are contained, the state |<i>^ can then be written 
as a product: 

1$) = |*ort.(l, 2, ..N)) \^sp^n{l, 2, ..N)) (51) 

with: 

\^orbXl.2,..N)) = \ui{l) U2{i)..UN{N))) (52) 

(particle 1 is inside measurement box 1, particle 2 inside measurement box 2, etc.) and: 

|^'sp„(l,2, ..iV)) = |1 : +; ...Af+ : +;Af+ + 1 : - N : -) + permutations (53) 

In this spin state, the spin up orientations and the spin down orientations are distributed in all 
possible ways among the numbered particles. 

If we wish, we can completely ignore the factorized orbital state and consider only the spins. For 
instance, with iV+ = 2 and = 1, the spin state reads as: 



1 

or, in a more condensed notation: 

1 

On the other hand, if = iV_ = 2, this spin state is, with the same notation 



^sp^nil, 2, 3)) = — [|1 : +; 2 : +; 3 : -) + |1 : +; 2 : -; 3 : +) + |1 : -; 2 : +; 3 : +)] (54) 



^sp^n{l,2,3)) [l-H, +,-) + ! + ,-, +) + !-,+,+)] (55) 



\^sp^nil,2,3A)) = 71 [! + ,+,-,-) + I + , -,+,-) + !-:+,-:+)+ ,5g. 
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etc. Such spin functions having equal amphtudes for all permutations belong to the category of W-states 
[5T1 [321 [33]. They can also be seen as the various M components of the "ferromagnetic states" [31] 
I J, M > of total angular momenturrF^ with the maximal possible value J = N/2. 

Our study shows that W-states are directly related to double Fock states and suggests a method to 
create them: start from a spin condensate and then perform a preliminary localization of the particles. 



7.2 Recovering the results obtained with identical particles 

We now check that a calculation with numbered spins allows us to recover our preceding results, without 
having to worry about orbital variables and symmetrization. For spin numbered j , the projector over the 
eigenstate corresponding to a result rjj = ±1 for a measurement along azimuthal direction (pj is: 

^,r"(^.) = \ [h + I + e^^^aj)] (57) 

with the usual notation cr^ for the angular momentum operator of spin j and Ij for the identity operator, 
which is the sum of the projectors over the two spin states 1+) and |— ): 

l. = l+>(+l, + l-X-l, (58) 

The probability of a sequence of a results 771, 772, ..t^at for measurements along polar directions Lpi, ip2, 
...(fiN is then proportional to: 

N 

{^sp^n{l,N)\l[-[\+){+\, + \-){-\, + |(e-'^^a+ + e'^^a-)] |vI/,,„(l,Ar)) (59) 

In this expression, each factor of the product of N brackets contains four terms, each with an operator 
that gives non zero only if it acts on one of the two states |+) and |— ); if a given choice among these 
four terms is made inside each bracket, a non-zero result is obtained for only one state for the N spins. 
For instance, if |— ) (— 1^, e~''^^cr^, e^'^^'a^ , and |-|-) (-I-I4, etc. are selected, the spin state has to be 
| — , — , +, +, ...). To obtain a non-zero result, a first condition is then that this state must have a non-zero 
component in the ket I'^spin)- To ensure that this condition is satisfied, we can multiply all |-|-) (-|-| 's and 
cr^'s by e*A, all |— ) (— |'s and cr+'s by e"'^^ and calculate the integral of the function F{A) obtained in 
this way by: 

'''^ ^^KN+-N-)Ap^^) (60) 

But a second condition is that the product with the bra (^spm(l, --i-^)! niust not vanish either, which 
is the case if the effect of the successive tT+ and cr^ operators flips the same number of spins in each 
direction. To ensure this, we multiply all (t+ 's by e*'*', all the (T~'s by e~'^^, and introduce a second 
integral over A: 

If these two conditions are satisfied, one always obtains a non-zero result, actually always the same number 
since all non-zero components of the state vector are equal. Finally, the probability is proportional tcF^: 



^''This property is only true for 1/2 spins; the proof can easily be obtained by recurrence, by applying several times 
operator J~ onto | J,J >. 

^"^Factors 2 disappear because, for instance, cr~ \ + >= 2 | — >. 
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or: 



/ — cos(iV+ - A^_)A / — [cosA + 77icos(A- A- (^j 
J-TT 27r 27r -^^x 



(63) 



which is identical to (fT^ with a trivial change of integration variable. We have therefore recovered the 
results obtained previously with identical particles, but with numbered spins, as usual in calculations of 
Bell inequalities violations. Our results are therefore not limited to Bose-Einstein condensates; we have 
a systematic way to go from identical to distinguishable particles. When the number of measurements 
M is less than the number of particles iV, a summation over the results r^j of the N — M unperformed 
measurements provides the probability, as in § 12.11 



7.3 Triplet state 

In the case N = M = 2, the initial state is the triplet state: 

1 r , 



I 4'spi„(l,2) >= 



71' 



and we obtain: 



or: 



— / — [cos A + ?7i cos (A' - [cos A + ?72 cos (A' - (^2) 



^ 1 + mm COS (<^i - (^2) 
or, after a normalization to 1 of the sum of the four different probabilities: 

1 



Pr. 



[1 +??l'72C0s((y9i -ip2)] 



This is the usual result, which can also be written as 

dX' 



P. 



2-K 



1 + 77i\/2cos (v3i - A') 1 + \/2?72 cos (i^2 - A') 



(64) 

(65) 
(66) 

(67) 
(68) 



In this expression, we see that the brackets inside the integral can indeed become negative, allowing a 
possible violation of the Bell inequalities. 



8 Conclusion 

Transverse spins measurements on double Fock states provide an interesting case where one can calculate 
exactly the predictions of quantum mechanics in all experimental situations, even if the measurements 
depend on many parameters. Another interesting flexibility arises from the choice of the two functions 
A and S, which can be defined in different ways; depending on this definition, the physical quantity on 
which locality is tested is microscopic, macroscopic, or intermediate; one can in this way study in detail 
the emergence of local classical properties of physical systems from microscopic non-locality, within the 
formalism of quantum mechanics, as a function of all the parameters of the experiment. 

One often rightly emphasizes that quantum entanglement is the essential ingredient of violations of 
local realism by quantum mechanics; the best known example of strongly entangled quantum states are 
the GHZ/NOON states, which indeed lead to strong violations. It is nevertheless interesting to note 
that this maximal entanglement is not a necessary condition; in fact, the minimum correlations that are 
compatible with Bose-Einstein statistics are already sufficient to lead to strong quantum non-local effects. 
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Conversely, this does not mean that statistical effects are a necessary condition for violations either; in 

fact, wc have shown that the same effects can be obtained with distinguishable spins in states belonging 
to the family of W states. A common property of all our results is that, in all cases, it is essential to 
perform the measurements on all particles; if a single one is missed, the violations disappear. When this 
condition is fulfilled, one reaches situations where the quasi-classical image of the Anderson phase is not 
always sufficient to reproduce the quantum results. 

Acknowledgments: Laboratoire Kastler Brossel is "UMR 8552 du CNRS, de I'ENS, et de I'Universite 
Pierre et Marie Curie" . 
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Appendix: angles of measurement 



To find the maximal violations of the Bell inequalities we used a numerical routine that produced the 
angles at which the violation occurred as well as the value of the violation. In the cases shown in Fig. [3] 
where A (and B, A' , and B' as well) is a product of results of measurements all at the same angle, we have 
noted previously that the angles at maximum occur in a fan arrangement with h — a = a — b' = a' — h = x, 
and a' — b' — 3% (in this appendix we simplify the angle notation by writing, for example (/?a = a). In 
Fig. m we give % for the same cases as treated in Fig. [3] For the case of P = N /2 the angle x drops off 
as 1/VN as shown analytically with the Gaussian approximation of Eq. 
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Figure 4: Angles of measurement corresponding to the violations shown in Fig. [21 The angle x is defines 
the fan described in the text. 

In the treatment of the "semi-mesoscopic" measurement described by Eq. (|151) , the angles also make 
a fan with x = ^^/^ for the one case = 6 that gives a violation. However, when we analyze the 
inequality described in Eq. (|44ll . the fan becomes distorted. We describe each of the situations of Table 
I individually. For = 4 all eight of the angles are distinct, but spread out in a pattern: Starting at a' 
and moving in order of a' ,d' ,b,c,a,d,b' ,c' we proceed alternately in steps of Ai and A2 where Ai — 0.458 
and A2 = 0.326 (with Ai+ A2 = 7i'/4). Thus the separation between any angle and its prime is 7r/2. 
The pattern is shown in Fig.[5l^a). 



For A^ = 6, we use one measurement each for A and for A' , but the product of two measurements 
for B and for B'; similarly for C and V and their primes. The two angles for B collapse to a single angle 
with the same holding for B', V, and T>' . We find a = c; if we take these as the origin, then we have a' at 
— 3A; b and d' at —A; b' and d at A; and c' at 3A where A — tt/8. The spread of the whole fan is then 
TT. See Fig. Mh) 

For A^ = 8 and 12, the angles collapse to just four distinct values. For A^ = 8 we have each letter A, 
B, A' , and B' being a product of two experimental results at the same angle. If a' = d' is the origin, then 
we move up in steps of Ai, A2, and Ai to b = c, a = d, and 6' = c', respectively, where Ai = 0.4533 and 
A2 = 0.1738. as seen in Fig.EJc). With A^ = 12, each letter represents the product of three experimental 
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Figure 5: Angles of measurement corresponding to the inequalities of Eqs. and Patterns (a) 

through (d) show angles for = 4, 6, 8, and 12, respectively, corresponding to the violations of Table I 
for inequality (|44p . Part (e) shows the angle pattern for the violation of Eq. ((45|) . 



results at the same angle. The arrangement is the same as for N — 8, except that the separations are 
reduced to Ai = 0.3741 and A2 — 0.0685, as shown in Fig. [^td). The Gaussian approximation shows the 
angle separations to be dropping as 1/ ^/N. 

Finally we have the inequality of Eq. (|45p . For N = 6, each letter represents just one measurement. 
The angles resulting in a maximum are all distinct but come in four evenly spaced groups of three as 
seen in Fig. O^e). If the set {a', e' , d'} is at {0, 0.159, 0.209}, then we move up by 7r/4 to the next triplet 
of {6, /, c}; 7r/4 to the next set of {a, e, d}; and another 7r/4 to {&', /', c'}. 
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